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ABSTRACT 

We explore the Gnedin-Ostriker suggestion that a post-Big-Bang 
photodissociation process may modify the primordial abundances of the light 
elements. We consider several specific models and discuss the general features 
that are necessary (but not necessarily sufficient) to make the model work. 
We find that with any significant processing, the final D and 3 He abundances, 
which are independent of their initial standard big bang nucleosynthesis 
(SBBN) values, rise quickly to a level several orders of magnitude above the 
observationally inferred primordial values. Solutions for specific models show 
that the only initial abundances that can be photoprocessed into agreement 
with observations are those that undergo virtually no processing and are already 
in agreement with observation. Thus it is unlikely that this model can work for 
any non-trivial case unless an artificial density and/or photon distribution is 
invoked. 

Subject headings: cosmology: theory - early universe - gamma rays: theory - 
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1. Introduction 



The standard model of big bang nucleosynthesis (SBBN) ( [Peebles 1966 ; Wagoner" 



Fowler, fc Hoyle 1967] ; [Schramm fc Wagoner 1977] ; |Boesgaard fc Steigman 1985| ; |Walker 



et al. 199l| ) makes well defined predictions of the primordial abundances of the light 



elements. These predictions agree reasonably well with the commonly accepted values of 
the primordial element abundances as inferred from observation ([Walker et al. 19911 ; [Bathes 



Boyd, fc Mathews 1993| ; [Smith, Kawano, fc Malaney 1993|) at a baryon-to-photon ratio of 



about i|~3x 10 -10 . However, this baryon density corresponds to only a few percent of the 
closure density of the universe. This result has motivated a number of attempts to reconcile 
the observed big bang abundances with a model in which rj is much larger. 

In this context, it was suggested by |Gnedin fc Ostriker 1992| (GO) that a population 
of black holes may have formed from an early generation of massive stars. The collapse 
of the stars to black holes does not significantly contaminate the primordial abundances. 
The black holes then accrete matter and produce a photon bath which further processes 
the remaining primordial material. GO chose initial SBBN abundances at two values of the 
baryon-to-photon ratio which are higher than that which agrees with any of the abundance 
determinations. Thus, they start with an excess of 4 He and 7 Li and a deficit of D and 3 He. 
The 4 He and 7 Li are then processed into the lighter elements as a function of the total 
photon energy input into the system. Although GO are not able to process any significant 
amounts of 4 He, they claim to be in agreement with the observed primordial abundances at 
the 3cr level. (See also a more recent paper QGnedin, Ostriker fc Rees 1995|) ). 



In this paper we investigate in detail general models of nucleosynthesis induced by 
post-Big-Bang photoprocessing and determine the conditions under which the observed 
light element abundances can be reproduced. As we will show, we are unable to find a 
non-trivial solution where the initial abundances do not already agree with observation. 
The reason for this is straightforward. The SBBN abundance of 4 He is larger, by some 4 
or more orders of magnitude, that those of D and/ or 3 He (especially at high values of if). 
It is therefore not possible to destroy enough 4 He to achieve consistency with observations 
without, at the same time, grossly overproducing pregalactic D and/or 3 He. Further, at 
high r], 7 Li is significantly overproduced and it is difficult to destroy enough 7 Li without 
destroying too much 4 He (and, consequently, overproducing D and 3 He). Our detailed 
calculations will provide support for this simple overview. 
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2. The Photon Spectrum and Energy- Weighted Cross Sections 



In the present scenario, the primordial elements are assumed to be bathed in a 
photon field for a given exposure time and photon intensity. In order to calculate the 
nucleosynthesis, it is necessary to know the energy spectrum of the radiation emitted by the 
black holes. Since such spectra are not observed directly, we look to other accreting objects, 
such as quasars, which may have similar spectra. Measurements have been made of NGC 
4151 ( |Perotti et al. 1981"! [Baity et al. 1984] ) and other quasars ( |von Ballmoos, Diehl, and 



Schonfelder 1987|) in energy ranges from 2 keV ( [Baity et al. 1984]) up to 100 MeV ( [Trombka 



et al. 1977| ; [Fichtel, Simpson, fc Thompson 1978| ). 



The photon number spectrum can be parameterized by two power laws, one describing 
the data below 2-3 MeV and with the form E~ 1S ( Baity et al. 1984] ), the other describing 



the data above 2-3 MeV with the form E~ 2 - 7 flRothschild et al. 1983| ). We are only 



concerned with energies above 2.2 MeV, the binding energy of the deuteron, and thus only 
use the higher energy power-law. It should be noted that |Gncdin fc Ostrikcr 1992| give a 



more complicated expression for the 7-ray spectrum in order to describe the full energy 
range. As we shall see, the choice of 7-ray spectrum does not change our conclusions. 

Cross sections for (7, X) reactions, where X = n, p, d, t, 3 He, or a have been measured 
over most of the energy range of interest. In order to describe the photoerosion process, we 
first calculate the energy-weighted cross sections (Poyd, Ferland fc Schramm 1989|) 



/ \ ^2MeV a i. E ) E 2 ' 7dE m 

\ a l = F-^dE ■ {) 

h.2MeV & a ^ 

Because of the power-law weighting factor, the energy-weighted cross section is dominated 
by contributions from the excitation function near threshold. Sparse data exist for 
photodissociation reactions on the lithium isotopes near threshold, contributing up to 
a 50% systematic uncertainty in the energy-weighted cross section. Cross sections for 
photodissociation reactions on 4 He, 3 He, and D are well determined from threshold up 
to several tens of MeV. The errors in these measurements are sufficiently small (< 20%) 
that our conclusions will not be affected by improved measurements of these cross sections. 
Table |2| contains a summary of the calculated energy-weighted cross sections. 
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3. Photoprocessing of the Elements 



We can describe the evolution of the element abundances locally by 

1 d^vi) = ^ ( ^ } _ (]Vi) (2) 

where Ni = Ni(f,t) is the number density (at location r and time t) of nuclide i and 
= 0(f) is the time-independent photon flux. We can further generalize Eqn. [2] by writing 
the number density as Ni(f,t) = fi (0(r)t) p(f), where p{r) is the local matter density. For 
a time-independent p(r), Eqn. ^| can be rewritten with the fi replacing the Ni and yielding 
solutions which are the same to within a dimensional constant. 

In general, Eqn. ^| might also include /3-decay terms which act as both production 
and destruction mechanisms. In our case the only unstable nuclides of importance are 3 H 
and neutrons. Reaction cross sections (and thus the energy-weighted cross sections) for 
photoproduction and photodissociation of 3 H are very similar to those for 3 He (see Table |2|). 
For neutrons the situation is even simpler since all are assumed to decay to protons without 
undergoing nuclear reactions. Therefore we treat 3 H as 3 He and neutrons as protons in our 
calculations. 

The simplest model is to adopt a uniform isotropic 7-ray distribution (i.e. is 
constant) and a uniform matter density distribution in Eqn. |2|. Such a distribution could 
be produced by a universe populated with a high density of black holes which tend to 
smear out the 1/r 2 dependence of the photon flux from an individual source. Eqn. |2] then 
describes a set of coupled differential equations which can be evolved forward in time as a 
function of the "exposure" (<f)t) for each element. Furthermore, the solutions obtained from 
this simple model can be modified to include effects from non-isotropic photon fluxes or 
non-uniform density distributions. These effects are folded into the solutions to Eqn. || by 
a weighting procedure of the form 



/ Ni \ J °° w(x)Ni(x)dx 



- 



(3) 



\N H / J ca w(x)N H (x)dx 

where x = (p(r)t is the position-dependent exposure. Note that with w(x) = 5(x — 4>(r)t) 
we recover the original local solutions. 

We examine the effects of several functional forms for w(x). In the first case we 
assume that there are two separately processed regions. Then w(x) is simply a sum of 
two S functions. In the second case, we have assumed a uniform population of black holes 
surrounded by a uniform density of matter. Processing of the matter is dominated by the 
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nearest black hole. Thus the final element abundances are given by a spatial averaging of 
the processed material. In the third case, we consider not only a distribution of black holes 
but a density distribution which varies as 1/r 2 so as to produce a flat rotation curve. From 
these examples, we reach some general conclusions about the necessary (and artificial) form 
w(x) must have in order to obtain a viable model. 

The criterion for these models to be viable is that the calculated abundances are 
within the limits determined by the observed light element abundances (see Table |1]). We 
have chosen realistic but generous (2a or greater) limits so that if a model does fail, the 
failure is due to the model not to the choice of limits. The standard model of big bang 
nucleosynthesis is in agreement with our limits for 2.5 x 10~ 10 < r] < 4.6 x 10~ 10 . 

In each of the scenarios mentioned above, we have presented the solutions in two 
different ways. First, we use initial abundances from standard big bang nucleosynthesis 
which allows us to parameterize our four initial abundances in terms of the baryon-to-photon 
ratio 7]. We can then follow the evolution of the abundances as a function of the exposure, 
<f)t. Second, since it can be argued that an assumption of the existence of primordial black 
holes negates the SBBN model and necessitates the use of the inhomogeneous model, 
we invert the solutions to Eqn. |2] for each of our scenarios. Then using the observed 
abundances, we can trace the evolution backwards in time to solve for the primordial values 
as a function of the exposure. Thus we avoid making any assumptions about the isotopic 
content of the pre-photoprocessed universe. 



4. Results 

4.1. Isotropic photon flux 

We have calculated the abundances of the light elements relative to hydrogen for 
isotropic exposures < (fit < 10 3 photons/10 -27 cm 2 starting with initial abundances 
determined from the SBBN model with 10 -10 < i] < 10~ 8 . Fig. [I] shows the evolution of the 
abundances as a function of the exposure, (fit. Horizontal lines represent the observationally 
inferred abundances. In Fig. |2] we map out the regions of (fit versus 77 space where the 
individual abundances agree with observation. At no values of rj and (fit can the processed 
abundances of 4 He, 3 He, and D be reconciled with the observed limits except in the region 
where the unprocessed abundances already essentially agree with observation. In fact, this 
scenario can be ruled out by just looking at the 4 He and 3 He abundances. Since most 
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4 He is being converted directly into 3 He and the 4 He abundance is so much greater than 
3 He, any small change in the 4 He abundance is accompanied by a very large change in the 
3 He abundance. This is clearly seen in the figures to be independent of the initial 3 He 
abundance (and therefore also independent of 77). Deuterium is constrained to rise more 
slowly than 3 He because its production is predominantly from 3 He through the 3 He(7,p)D 
reaction. Therefore, deuterium production cannot proceed until a significant quantity of 
3 He has been created. 

At large exposures, when most of the 4 He has been destroyed, the 3 He abundance 
begins to decrease, eventually reaching agreement with observation. However, this large <j)t 
scenario is ruled out because the 4 He and D abundances are well outside the observational 
bounds. 

We have inverted the solutions to Eqn. ||] in order to calculate the possible combinations 
of initial abundances for which this model will work. We are therefore no longer constrained 
by SBBN. Solutions for the initial abundances, which when processed for a given <pt yield 
light element abundances in agreement with observations, are found only for <pt < 9.1 x 10~ 3 
photons/1 0~ 27 cm 2 . Therefore, no significant processing occurs. Indeed, the initial 4 He 
and 3 He abundances are constrained to be within the observational limits. The initial D 
abundance can be less than the observed lower limit by only ~ 5% and cannot be more 
than the observed upper limit. The initial 7 Li abundance can be no more than 0.3% above 
the upper limit. 



4.2. Two-zone model 

We have investigated the possibility of reconciling the helium and deuterium isotope 
abundances by assuming a two-zone model, i. e. w(x) is a sum of two 5 functions. This 
provides a simplified description of a mixture of material processed through different 
exposures. In this case, space is divided into two regions, one of which undergoes very little 
or no processing and the other undergoes a great deal of processing. In the former case, 
deuterium and 3 He agree with observation and 4 He is overabundant. It is also necessary 
that rj < 10 x 10~ 10 so that deuterium is not underproduced (see Fig. [I]). The latter zone 
undergoes significantly more processing such that most of the 4 He is destroyed. Also, the 
3 He and deuterium which have been created from the processed 4 He have been destroyed to 
such an extent that they once again are close to the observational limits. Then by adjusting 
the relative volumes of the two zones (such that the larger volume is only slightly processed 
and the smaller volume is highly processed), one can force agreement between the processed 
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4 He and the observational limits. In the larger mostly-unprocessed zone the 3 He abundance 
constrains <pt < 10~ 2 . The 4 He and 7 Li are therefore essentially unprocessed. However, 
since 7 Li is overproduced by a greater factor than 4 He, and the zone-mixing reduces the 7 Li 
by the same fraction as 4 He, the 7 Li is still overproduced in this model. The only way to 
reconcile the model with observation is to start with primordial abundances for 7 Li, 3 He, 
and deuterium which agree with the observational limits already. Even then, the primordial 
4 He cannot be more than a few percent overabundant due to the constraints placed on rj by 
7 Li. 



4.3. Uniform black hole distribution 



Although the two-zone model shows clearly why agreement cannot be attained between 
the primordial abundances and observation, it is too simple a model from which to infer 
general conclusions. In order to investigate a more physical picture, we have studied a 
model that has a uniform distribution of black holes that are assumed to be weak enough 
and/or far enough apart that processing of the primordial material is dominated by the 
closest black hole. In order to account for the spatial variation of the photon flux, we 
modify Eqn. ^| by setting <f>(r) = L^/Airr 2 , where L 7 is the total 7-ray luminosity from 
one accreting black hole. We can then write the spatially-averaged abundance relative to 
hydrogen as 



Nh 1 It:: N H r 2 dr flte f H p(r)r 2 dr ' 

We examine two cases in detail. In the first case, we assume that p(r) is constant, i.e. 
the matter is distributed homogeneously. Then Eqn. § can be written more conveniently in 
terms of one parameter x = to yield 

Ni \ [* max fi (x) x~ 5/2 dx , 



N H I f* max Ih (x) x~ 5 / 2 dx ' 

Alternatively, if we assume that the density of matter falls off with distance from the 
central source as p(r) = Po(f) 2 (which will yield the flat rotation curves seen in galaxies) 
then Eqn. f| becomes 

/ Ni\ J* max fi (x) x~ 3/2 dx , , 



- 



\ Njj I r max f H (x) x~ 3 / 2 dx ' 
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Figures |3| and |] show the dependence of the abundances on x min (with x max = oo so 
that R m i n = 0) for three values of r\. Figures || and ^| illustrate the regions of r\ — x min space 
where the spatially-averaged abundances agree with observation. It is clearly seen that there 
is no region of parameter space where the models agree with the observed abundances of 
4 He, 3 He, and D simultaneously except where the standard model (for low baryon density) 
already works. We have also searched for solutions for x min < x max < oo, where processed 
material with exposure greater than x max is assumed to be swallowed by the black hole and 
thus does not contribute to the spatial average. No non-trivial solutions were found. 

As with the isotropic photon flux scenario, we have inverted the solutions to Eqns. |5| 
and ^| in order to calculate the possible combinations of initial abundances for which these 
models will succeed. Solutions for the initial abundances, which when processed yield light 
element abundances in agreement with observations, are found only for <pt < 3 x 10~ 3 and 
cf)t < 8 x 10~ 6 photons/10 -27 cm 2 for the uniform density and the flat rotation curve models, 
respectively. Again, no significant processing occurs and the initial 4 He and 3 He abundances 
are constrained to be within the observed limits. The initial D abundance may be less than 
the observed lower limit by only ~ 5% and it cannot be higher than the observed upper 
limit. The initial 7 Li abundance can be no more than 0.3% above the upper limit. 

It might be thought that the detailed shape of the photon spectrum will affect the 
present results. However, changing the power law index or using the parameterization of 
GO causes only small changes in the energy-weighted cross sections. In the GO case, the 
energy-weighted cross sections are increased by no more than a factor of 4 even when using 
the most extreme distortion of the spectrum. Furthermore, changes in the spectra change 
all the energy-weighted cross sections in the same direction, so the effect on resulting 
abundances is less than that on the cross sections. Since our conclusions are based on 
qualitative conflicts between observations and the model predictions, it is not likely that 
uncertainties in the photon spectrum will affect our conclusions. 



4.4. A general solution? 

So far we have provided examples of models that don't work. Is it possible to construct 
a non-trivial model that will work? That is, can we construct a weighting function w(x) 
(see Eqn. |3|) that will reconcile the solutions to Eqn. |2| (shown in Fig. HD with the observed 
limits? In constructing such a model, two important features must be kept in mind. First, 
the weighting function must preferentially mix regions such that a larger fraction of 7 Li is 
destroyed than of 4 He. Since the destruction cross section for 7 Li is higher than for 4 He, 



- 9- 



this may be possible. Second, the weighting function must avoid appreciable contributions 
from the region 10~ 2 < <pt < 5 x 10 2 in order to avoid overproducing 3 He and D. As a result, 
such a function is forced to be of the form shown in Fig. [T[ This bimodal distribution 
is decidedly ad hoc! But, even with a weighting function of this type, it is difficult (and 
perhaps impossible) to reproduce the observed primordial abundances of D, 3 He, 4 He, and 
7 Li. 



5. Conclusions 

We have investigated several scenarios in which the primordial abundances of the 
light elements undergo photoprocessing due to a population of 7-emitting objects such as 
accreting black holes (GO). In all cases, the primordial 4 He abundance is constrained to 
agree with observation since even a small amount of processing of 4 He increases the 3 He 
and D abundances by several orders of magnitude. At large exposures, when destruction 
of 3 He and D has occurred, their abundances relative to hydrogen become independent of 
initial values and cannot simultaneously be reconciled with observation. A general solution 
which would correct for this effect must preferentially weight the processed regions in such 
a way as to be pathological. 

In short, we find no non-trivial solution to any of the several photodissociation models 
considered for either SBBN or inhomogeneous BBN abundances corresponding to large 
values of rj. Post BBN photoprocessing of the light elements cannot weaken the upper 
bound on the universal density of baryons. 

This work was supported in part by National Science Foundation grant PHY92-21669 
and Department of Energy grant DE-AC02-76ER01545. 
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Fig. 1. — Abundances relative to hydrogen of the light elements as functions of exposure, 
4>t. The initial abundances are calculated from SBBN with 7710 = 3 (solid curve), 7710 = 10 
(dashed curve), and r] 10 = 100 (dashed-dotted curve). Horizontal lines indicate observational 
limits on the primordial abundances. 

Fig. 2. — Parameter space for the isotropic photon flux model. Hatched areas indicate 
agreement between the calculated and observed elemental abundance. Vertical hatching is 
D, hatching with a positive slope is 3 He, and hatching with a negative slope is 4 He. The 
region where 7 Li agrees with observation is not shown, however it also overlaps the region of 
agreement between the other three light elements. 

Fig. 3. — Abundances relative to hydrogen of the light elements when spatially averaged 
over a uniform density distribution and a photon flux which decreases as 1/r 2 . The abscissa 
is x min , where x = Lt/4nR 2 . The initial abundances are calculated from SBBN with r] w = 3 
(solid curve), i] W = 10 (dashed curve), and r] W = 100 (dashed-dotted curve). Horizontal 
lines indicate observational limits on the primordial abundances. 

Fig. 4. — Abundances relative to hydrogen of the light elements when spatially averaged 
over a density distribution and a photon flux both of which decrease as 1/r 2 . The abscissa is 
Xmin, where x = Lt/AnR 2 . The initial abundances are calculated from SBBN with r] 10 = 3 
(solid curve), 7710 = 10 (dashed curve), and r7 10 = 100 (dashed-dotted curve). Horizontal 
lines indicate observational limits on the primordial abundances. 

Fig. 5. — Parameter space for the spatially- averaged model with a uniform density 
distribution. Hatched regions indicate where the calculated abundances of D, 3 He, and 
4 He agree with observation. Vertical hatching is D, hatching with a positive slope is 3 He, 
and hatching with a negative slope is 4 He. The only agreement for all three abundances 
occurs in the region of parameter space where the standard (low baryon density) model 
already works. The region where 7 Li agrees with observation is not shown, however it also 
overlaps the region of agreement between the other three light elements. 

Fig. 6. — Parameter space for the spatially- averaged model with a matter density which 
varies as 1/r 2 . Hatched regions indicate where the calculated abundances of D, 3 He, and 
4 He agree with observation. Vertical hatching is D, hatching with a positive slope is 3 He, and 
hatching with a negative slope is 4 He. The only agreement for all three abundances occurs 
in the region of parameter space where the standard (low baryon density) model already 
works. The region where 7 Li agrees with observation is not shown, however it also overlaps 
the region of agreement between the other three light elements. 
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Fig. 7. — Possible weighting function w(x) (thick line) superimposed over the abundance vs. 
exposure curve for 3 He. 



Table 1: Observational limits on primordial abundances. 
Abundance Limit Reference 
7 Li/H < 7 x 1(T 10 1, 2, 3 

Y p 0.215 - 0.244 1, 2, 4, 5 

3 He/H < 2 x 10~ 5 1 
D/H 1.2 - 10 xl0~ 5 1, 2 



References. — 1) Walker ct al. 1991 , 2) smith, Kawano, fc Malaney 1993 , 3) Pinsonncault, Deliyannis, fe 



Demarque 1992| , 4) |Balbes, Boyd, fe Mathews 1993| 5) plive fc Steigman 199E| 
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Table 2: Energy-weighted cross sections 

Reaction < o > (mb) Reference 

7 Li( 7 , n) 6 Li 0.0486 1 

7 Li( 7 , p) 6 He 0.0480 2, 3 

7 Li( 7 , d) 5 He 0.0118 2 

7 Li( 7 , t) 4 He 0.1760 2, 4, 5, 6 

6 Li( 7 , n) 5 Li 0.1748 1 

6 Li( 7 , p) 5 He 0.4334 2 

6 Li( 7 , d) 4 He 0.0047 7 

6 Li(7, t) 3 He 0.0147 2 

4 He(7, n) 3 He 0.0153 8, 9, 10 

4 He(7, p) 3 H 0.0167 8, 11 

4 He( 7 , d)D 7.7 x lO" 5 12 

4 He(7, pn)D 0.0013 13 

3 He(7, n)pp 0.0547 14 

3 He(7, p)D 0.0913 15 

3 H(7, n)D 0.0547 14 

3 H( 7 , p)nn 0.0788 14 

D( 7 , p)n 1.7473 16 



References. — 1) Dietrich & Berman 1988, 2) Junhgans et al. 1979, 3) Gregory, Sherwood, & Titterton 



1962|, 4) [Griffiths ct al. 196l|, 5) |Burzynski ct al. 1987j, 6) |5kopik et al. 1979|, 7 ) |Robcrtson ct al. 1984 8) 
Feldman ct al. 1990| , 9) |Ward et al. 1981] , 10) [Berman 1980|, 11) |Bernabei 1988| 12) |Barnes et al. 1987| , 13) 
Balestra et al. 1979| , 14) |Faul et al. 198l| , 15) fTiccioni ct al. 1973| 16) |Segre 1964 



D/H 



Y, 



o 



o 



o 



I 



I 



o 



IT" 

o 

o 

CO 



o 



o - 



o 



OJ 



1 1 






- 


- 1 

1 1 

- \ 1 




- 


\ 


\ 

\ 






\\ 


- 






\ -_ 

r 


1 1 


1 


, I , | i 



; He/H 





'Li/H 




